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1.  Introduction 


The  crystal  structure  of  silicon  carbide  (SiC)  is  dictated  by  its  polytype.  A  SiC  polytype  is 
represented  by  the  number  of  hexagonal  planar  SiC  double  layers  in  its  unit  cell,  with  an 
appended  C,  H,  or  R  denoting  a  cubic,  hexagonal,  or  rhombohedral  structure.  Over  200 
structural  polytypes  of  SiC  single  crystals  are  known  to  exist  (Bernstein  et  ah,  2005).  Cubic 
polytypes  are  labeled  P-SiC,  while  hexagonal  and  rhombohedral  polytypes  are  labeled  a-SiC. 

The  common  cubic  (P)  polytype  is  3C-SiC,  exhibiting  the  sphalerite  or  zinc  blende  crystal 
structure.  Technologically  relevant  polytypes  of  hexagonal  (a)  SiC  include  21 1  (wurtzite  crystal 
structure),  AH,  and  6H,  while  common  rhombohedral  polytypes  include  9 R,  12 R,  15 R,  and  21 R. 
Hexagonal  polytype  6//-SiC  is  the  focus  of  the  present  investigation;  crystals  of  6//-SiC  are 
often  referred  to  as  moissanite.  Polycrystalline  structural  ceramics,  such  as  SiC-B  and  SiC-N  are 
thought  to  consist  primarily  of  6//-SiC,  though  other  polytypes  may  be  present  (Shih  et  ah, 

2000).  Polytype  6//-SiC  is  non-centrosymmetric,  piezoelectric,  and  spontaneously  polar,  i.e. 
pyroelectric  (Qteish  et  ah,  1993).  Phase  transfonnations  at  ambient  temperatures  are  thought  to 
occur  at  very  high  pressures  on  the  order  of  100  GPa  (Vogler  et  al.,  2006)  though  fractional 
increases  of  6//-SiC  grains  relative  to  other  polytypes  have  been  observed  to  occur  at  somewhat 
lower  compressive  stresses  ( — 19—32  GPa)  in  polycrystalline  SiC  subject  to  ballistic  impact  (Shih 
et  al.,  2000). 

Defects  can  strongly  affect  mechanical,  thermal,  and  electrical  behaviors  of  ceramics,  such  as 
SiC.  Basal  slip  of  edge  dislocations  is  thought  to  be  the  dominant  mechanism  of  plastic 
deformation  in  hexagonal  SiC  polytypes,  with  dislocation  mobility  dependent  on  temperature 
and  electric  current  (Galeckas  et  al.,  2002).  Shockley  partial  dislocations  bordering  intrinsic 
stacking  faults  are  thought  to  be  most  relevant  in  hexagonal  SiC  (Samant  et  al.,  1998;  Blumeneau 
et  al.,  2003),  with  the  mobility  of  leading  partials  exceeding  that  of  trailing  partials  at  low 
temperatures  (Demenet  et  al.,  2000)  supporting  a  tendency  for  plasticity  to  be  accompanied  by 
generation  of  large  numbers  of  stacking  faults.  Twinning,  which  is  observed  in  p-SiC  (Van 
Torne,  1966)  does  not  appear  to  be  an  important  deformation  mechanism  in  hexagonal  polytypes 
except  under  high  temperatures  and  pressures  wherein  phase  transformations  occur.  Point 
defects,  such  as  Si  and  C  vacancies  (Shih  et  al.,  2000;  Bernstein  et  al.,  2005)  and  interstitials, 
such  as  those  generated  by  irradiation  (Gao  and  Weber,  2005)  may  affect  physical  properties. 
Primary  cleavage  planes  are  of  basal  and  prismatic  orientations  with  single  crystals  of  a-SiC 
apparently  less  resistant,  in  some  cases,  to  fracture  during  indentation  experiments  than  their 
polycrystalline  counterparts  (Kitahara  et  al.,  2001)  suggesting  that  cleavage  of  single  crystals,  in 
addition  to  grain  boundary  fracture,  is  important  in  failure  assessment  of  polycrystals. 

This  report  addresses  the  behavior  of  single  crystals  of  6/Z-SiC.  A  model  incorporating  coupled 
non-linear  thermo-electro-mechanical  behavior  is  developed  in  section  2  for  proper  interpretation 
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of  results  of  shock  physics  experiments  on  anisotropic  piezoelectric  single  crystals  (Graham, 
1972;  Clayton,  2010).  The  primary  content  of  this  report  consists,  in  section  3,  of  a  number  of 
tables  of  structural,  mechanical,  thermal,  and  electrical  properties — either  obtained  from  the  open 
literature  or  newly  computed — of  hexagonal  polytypes  of  SiC  single  crystals.  Compression  data 
for  polycrystalline  SiC  has  been  reviewed  extensively  elsewhere  (Dandekar,  2002). 

The  notation  of  non-linear  continuum  mechanics  is  used.  The  index  notation  follows  the 
Einstein  summation  convention.  Current  configuration  indices  are  written  in  lower  case  Latin, 
reference  configuration  indices  in  upper  case  Latin,  and  intermediate  configuration  indices  in 
Greek.  Superposed  -1  ,T,  and  •  denote  the  inverse,  transpose,  and  material  time  derivative, 
respectively.  Subscripted  commas  and  semicolons  denote  partial  and  covariant  differentiation, 
respectively. 


2.  General  Model  for  Elastic  Dielectrics  with  Defects 


The  framework  of  section  2  combines  treatments  of  crystal  plasticity  of  ceramic  single  crystals 
(Clayton,  2009a)  and  dielectrics  with  point  defects  (Clayton,  2009b)  and  may  be  applied  to 
describe  any  single  crystalline  dielectric  or  semiconducting  solid  containing  dislocations, 
stacking  faults,  and  point  defects.  Let  xa  =  x“(XA,t)  denote  spatial  coordinates  that  depend  on 
reference  coordinates  XA  of  a  material  particle  and  time  t.  The  deformation  gradient  in 
equations  1-3  is 


/7"  =xa  =  FEa  Fa  =  FLa  FPa 

r.A  Xa  r  ,ar.A  r  ,ar  .A  ’ 


j^La  pEa  j  1/3 


} 7a  -  jUlpPa 
1  A  J  1  .A  > 


(1) 

(2) 

(3) 


where 

FEaa  :  reversible  thennoelastic  deformation 
FPaA  :  isochoric  deformation  from  dislocation  glide 

J  :  volume  changes  due  to  defects  (Eshelby,  1954;  Clayton,  2009a) 
FLaa  :  total  lattice  deformation 

F“  :  total  residual  deformation 
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In  equation  4,  spatial  velocity  gradient  Lab  is 


l\ = v.:  =  pEaFE~ia, + FEar:FE-'{ +~jj 


.b  1  ^  .a^.p 

r 


=  phaF 


Ea  T-iE-\a 


+  F 


Ea 


Z*  i — ia  — i 

ys  mp 


V  i 


F 


E-la 


(4) 


where 

va  =  xa  :  spatial  velocity  or  particle  velocity 

Lpp  :  plastic  velocity  gradient  referred  to  the  intennediate  configuration 

y‘  :  slip  rate  defined  for  each  slip  system  i 

J‘a  :  unit  slip  direction  defined  for  each  slip  system  i 

wip  :  unit  slip  plane  nonnal  defined  for  each  slip  system  i 

gah  =  g„  »gA  :  spatial  metric  tensor 

The  metric  tensor  in  the  intermediate  configuration  is  assumed  Cartesian  for  simplicity,  i.e., 

Sap  ^ ap  ' 


In  equations  5-8,  Maxwell’s  equations  in  the  quasi-electrostatic  approximation  (Maugin,  1988) 
for  non-magnetic  solids  are  written  as 

sabcec.,b=  0, 


dda 


J 


dt 
d-a=P 


K  =  o, 


(5) 

(6) 

(7) 

(8) 


where 

ea  =  -(j)a  =  ~<t>.a  :  spatial  electric  field 

da  =  s0ea  +  pa  :  electric  displacement 
s0  :  vacuum  pennittivity 

p“  :  electric  polarization  per  unit  spatial  volume 

/"  =/“-(!/ p0)£abcbc.b  :  effective  electric  (displacement)  current  density 
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i a  :  free  electric  current  (vanishes  in  a  perfect  dielectric) 
juQ :  vacuum  penneability 

p  :  free  charge  density  per  unit  spatial  volume  (vanishes  in  a  perfect  dielectric) 
ba  :  magnetic  flux  density 

sabc  =  I” ^det(gaA)  eabc  :  =  alternator  tensor  or  permutation  tensor 


e  :  internal  energy  per  unit  mass 
qa  :  heat  flux  (i.e.,  conduction) 
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r  :  heat  source  per  unit  mass  (i.e.,  radiation) 

\//  =e- 9rj  :  Helmholtz  free  energy  per  unit  mass 

9  :  absolute  temperature 
77  :  entropy  per  unit  mass 

In  equation  15,  the  free  energy  per  unit  mass  is  of  the  functional  form 

Y  =  Y(EEaP,pa,9,Z,S,9),  (15) 


where,  in  equations  16-17,  the  elastic  strain  tensor  and  intermediate  polarization  vector  are, 
respectively, 


1 


1 


Kp  =  2  (Q  -<V> =^F  agabF  e  -<v> = 


—  rpEa  * 

Pa=F  .aPa- 


(16) 


(17) 


Also,  E, ,  g ,  and  9  are  the  scalar  dislocation  length  per  unit  volume,  stacking  fault  area  per  unit 
volume,  and  point  defect  concentration,  respectively.  In  equations  18-21,  constitutive 
relationships  for  stress,  electric  field,  entropy,  and  heat  flux  consistent  with  equations  13  and  14 
are 


-ap  _  jE  t-<  E-la  ab  E-\p  _  -  <^Y  ,  TErE-laX-  rE-\p5- 


Zap  =  J*F^ayt>F*-W  =  p  __L__  +  Jh Cb~laX p  C 


BEe 

UnaP 


-s  ’ 


7  = 


dpa  8pa 

dy 


where 


39 

r =~KaP  9. 

p  =  pJE  =(p/  9)sahc£aPx FEaaFEb/3FEcx  :  intermediate  mass  density 

es  =  FEasea  =  ~FEas(j).a  :  intennediate  electric  field 
-a  _  jEpE-ia^a  .  intermeciiatc  heat  flux 


(18) 

(19) 

(20) 
(21) 


9p  =  9.aFEap  ■  anholonomic  temperature  gradient 
K"ft  :  thermal  conductivity 
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3.  6//-SiC:  Structure  and  Physical  Properties 


In  equation  22,  the  Helmholtz  free  energy  is  written  as  follows  on  a  per  unit  reference  volume 
basis,  where  8()  is  a  reference  temperature: 

PoV  =  \c*sElf  E‘s  +i 

-P^E^e-e,)-  P„cre  log  2  iV  (#.  ? .  •») , 


where 


p0  =  p.J  =  pJEJ  =  (p/ 6)sabc£ABC  Fa4FhgFcc  :  reference  mass  density 


Equations  23-31  apply  for  isothermal  permittivity  (i.e.,  dimensionless  dielectric  constants)  saf, , 
standard  piezoelectric  constants  efx ,  pressure  derivatives  of  dielectric  constants  at  the  reference 
state  os"11  /  dp  ,  second-order  elastic  constants  at  fixed  electric  field  ,  thermal  expansion 
coefficients  at  constant  polarization  ai/f! ,  isobaric  specific  heat  per  unit  mass  cp  ,  Gruneisen 
numbers  raP  and  r  ,  and  isothermal  bulk  modulus  B  (Polian  et  ah,  1982): 

£aP  =  +  £0  '  Ajff  ’ 

ef=-A~lsAspx,  (24) 


dp  3Bs0 


(YxS; A-XaxA-xps  -  s08ap  +  A~lafi) , 


QuPxS  _  QaPxS  _  ^eaP ^>xS 

pap  _  CaPx5a  s. 


c=cv+(ap)2  — , 


pap  —  _ 
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(31) 


B  = 


C33(C11  +  C12)-2(C13)2 
Cu+C12-4C13+2C33 


Only  two  of  the  six  possible  coefficients  of  TaPxS  are  known  from  equation  25;  it  is  assumed 
here  for  simplicity  that  Ta,lz<i  =  TaP 8xd  .  Pyroelectricity  is  omitted  in  equation  22.  Differences 
between  isothermal  and  isentropic  elastic  coefficients  are  negligible  in  SiC  near  room 
temperature. 


In  equations  32-34,  relevant  dislocations  in  a-SiC  include  partial  dislocations  on  the  basal  plane 
formed  from  the  dissociation  reaction  J^[1210]  — »  /3[l  100]  +  X[0  1 10]  (Maeda  et  ah,  1988; 
Samant  et  ah,  1998;  Demenet  et  al.,  2000;  Blumeneau  et  ah,  2003).  Dislocation  kinetics  are 
dictated  by  (Demenet  et  al.,  2000): 


f=A 


exp 


f zQ} 


sgn 


m  = 


d  log  t  1 
d  logy' 


fr_] 

(32) 

(33) 

(34) 

where 


A  :  a  proportionality  constant  with  dimensions  of  strain  rate 

Q  :  activation  energy  that  may  depend  on  electric  current  (Galeckas  et  al.,  2002) 

m  :  strain  rate  sensitivity 

t‘c  :  intrinsic  slip  resistance  or  shear  strength 


In  equation  35,  residual  elastic  volume  changes  from  screw  dislocations  are  estimated  as 


AV/V  =  J -lx  — 
M 


dju  n 
dp  B 


(35) 


where 

p  :  effective  elastic  shear  modulus 

Ee  :  elastic  energy  per  unit  dislocation  line  length 


Analogous  fonnulae  for  edge  dislocations  and  point  defects  are  available  in  the  literature 
(Eshelby,  1954;  Holder  and  Granato,  1969;  Clayton,  2009a;  Clayton,  2010). 

Table  1  lists  parameters  of  the  6//-SiC  crystal  structure.  The  stacking  sequence  of  Si  or  C  layers 
in  the  [0001]  direction  (i.e.,  along  the  c-axis)  is  ABCACBABCACB...,  meaning  six  hexagonal 
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planar  layers  each  of  Si  and  C  comprise  a  single  conventional  unit  cell.  Tables  2  and  3  list 
fundamental  and  derived  thermal,  mechanical,  and  electrical  properties  of  a-SiC.  Material 
property  tensors  are  expressed  using  Voigt’s  condensed  notation.  Tables  4-7  list  properties  and 
supporting  references  regarding  defects  and  fracture.  Predictions  for  volume  changes  associated 
with  defects  (see  equation  35  and  other  equations  in  [Clayton,  2010])  are  shown  in  figure  1. 


Table  1.  Structural  features  of  6 //- S i C  (room  temperature  and  atmospheric  pressure). 


Parameter 

Value 

Remarks 

Reference 

a 

0.308  nm 

Lattice  parameter  (C-axis) 

Stockmeier  et  al.  (2009) 

c 

1.512  nm 

Lattice  parameter  (V-axis) 

fio 

0.0103  nm3 

Atomic  volume  </(,  =  V3a2c/24 

Point  group 

6  mm 

Dihexagonal  pyramidal  (HI) 

Arlt  and  Schodder  (1965) 

Space  group 

C^y4  (C6mc) 

Wyckoff  (1963) 

Table  2.  Thermophysical  properties  (room  temperature  and  atmospheric  pressure). 


Parameter 

Value 

Remarks 

Reference 

Po 

3210  kg/m3 

Mass  density 

Slack  (1964) 

CP 

652  J/kgK 

Specific  heat  at  constant  pressure 

Muller  et  al.  (1998) 

an 

3. 4x1  O'6  /K 

Thermal  expansion  coefficients 

Li  and  Bradt  (1987) 

a33 

3. 2x1  O'6  /K 

Pn 

2.3  MPa/K 

Thermal  stress  coefficients 

equation  27 

Pi  3 

2.1  MPa/K 

Cl 

1.10 

Gmneisen  coefficients 

equation  29 

r33 

1.00 

r 

1.06 

Gruneisen  parameter 

equation  30 

Ku 

390  W/mK 

Thermal  conductivity 

Burgemeister  et  al.  (1979) 

K3 3 

370  W/mK 

©D 

1200  K 

Debye  temperature 

Slack  (1964) 

©M 

3375  K 

Melting  temperature 

Merala  et  al.  (1988) 
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Table  3.  Elastic,  piezoelectric,  and  dielectric  properties  (room  temperature  and  atmospheric  pressure). 


Parameter 

Value 

Remarks 

Reference 

c„ 

501  GPa 

Second-order  elastic  constants; 
Brillouin  scattering; 

6/f-SiC 

Kamitani  et  al.  (1997) 

Cl  2 

112  GPa 

C\  3 

52  GPa 

C33 

549  GPa 

C44 

161  GPa 

dCn/dp 

3.8 

Pressure  derivatives  of  second- 
order  elastic  coefficients; 
atomic  model; 

2//-SiC 

Davydov  (2004) 

dCu/dp 

4.0 

dC  1 3/ dp 

4.0 

dCu/dp 

3.8 

SC44/  dp 

-0.2 

t* 

194  GPa 

Isotropic  elastic  constants; 

Voigt  averages;  6 / /- S  i C 

Blumeneau  et  al.  (2003) 

V 

0.161 

B 

222  GPa 

equation  3 1 

B 

230  GPa 

Static  compression;  polycrystal 

Bassett  et  al.  (1993) 

dB/dp  =  B' 

4.0 

Static  compression;  polycrystal 

Bassett  et  al.  (1993) 

3. 1+0.3 

Shock  compression;  polycrystal 

Sekine  and  Kobayashi 
(1997) 

3.9 

Isotropic  elasticity,  dCa/Jdp 

2.5 

dB/dp  =  2.T+1/3 

Slater  (1939) 

dju/dp  =  fu' 

3.4 

Isotropic  elasticity  at  fixed  v 

e\5 

-0.198  C/m2 

Piezoelectric  constants; 
atomic  model;  6//-S1C 

Mirgorodsky  et  al. 

(1995) 

e3i 

-0.200  C/m2 

e33 

0.398  C/m2 

£ll 

9.66 

Static  dielectric  permittivity; 
experiment;  6//-SiC 

Patrick  and  Choyke 
(1970) 

£33 

10.03 

0£\ \!dp 

-0.0218 /GPa 

Pressure  derivatives  of 
permittivity; 

DFT  calculation;  4//-SiC 

Karch  et  al.  (1996) 

dsa/dp 

-0.0232  /GPa 

An 

13.0  m/nF 

Inverse  dielectric  susceptibilities 

equation  23 

12.5  m/nF 

a15 

2.58xl09  C/(Fm) 

Stress-polarization  coefficients 

equation  24 

An 

2.50xl09  C/(Fm) 

A3  3 

-4.98xl09  C/(Fm) 

r,  1 

-5.67  m/nF 

Piezo-optic  coefficients 
(electrostriction) 

equation  25 

^33 

-5.49  m/nF 
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Table  4.  Slip  systems  and  Burgers  vectors. 


Type 

Direction 

Plane 

Burgers  b 

1*1 

Planar  spacing  d 

Basal,  frill 

<  1210  > 

(0001) 

1/3[1210] 

0.308  nm 

0.252  nm  (c/6) 

Basal,  partial 

<iToo> 

(0001) 

1/ 3[1  Too] 

0.178  nm 

0.252  nm  (c/6) 

Prism,  full 

< 1210  > 

{io  To} 

1/3[1210] 

0.308  nm 

0.267nm(aV3/2) 

Table  5.  Fracture  properties. 


Parameter 

Value 

Plane 

Remarks 

Reference 

Yf 

8.6+0.4  J/m2 

(0001) 

Cohesive  energy; 
indentation  experiment 

Kitahara  et  al.  (2001) 

5.6+0.5  J/m2 

{ioTo} 

4.0+0.2  J/m2 

{T2T0} 

37.7  GPa 

(0001) 

Theoretical  strength 
/ub 

X  p  — 

2nd 

Frenkel  (1926) 

35.6  GPa 

{10T0} 

Table  6.  Dislocation  properties  (unless  noted  otherwise:  basal  dislocations,  -room  temperature). 


Parameter 

Value 

Defect  type 

Remarks 

Reference 

Ee 

21.1  nJ/m 

Full  60° 

Elastic  energy;  isotropic 

Blumeneau  et  al.  (2003) 

6.1  nJ/m 

Partial  30° 

Elastic  energy;  anisotropic 

7.2  nJ/m 

Partial  90° 

Elastic  energy;  anisotropic 

18.5  nJ/m 

Edge  loop 

Elastic  energy;  isotropic 

Ec 

0.72  nJ/m 

Partial  30°  (Si) 

Core  energy;  DFT;  2//-SiC 

Blumeneau  et  al.  (2003) 

0.74  nJ/m 

Partial  90°  (Si) 

Core  energy;  DFT;  2//-SiC 

1.25  nJ/m 

Partial  30°  (C) 

Core  energy;  DFT;  2//-SiC 

1.81  nJ/m 

Partial  90°  (C) 

Core  energy;  DFT;  2//-SiC 

WSF 

2. 5+0. 9  mJ/m2 

(0001)  plane 

Stacking  fault  energy;  -1800K 

Maeda  et  al.  (1988) 

4.5  GPa 

Unknown 

Shock;  polycrystal  @  HEL 

Feng  et  al.  (1998) 

3. 3+1.0  GPa 

Unknown 

Shock;  polycrystal  above  HEL 

Vogler  et  al.  (2006) 

6.8  GPa 

Unknown 

Static  compression;  polycrystal 

Zhang  et  al.  (2002) 

4. 5+0. 4  GPa 

Unknown 

Indentation;  single  crystal 

Shim  et  al.  (2008) 

0.45  GPa 

Partial  basal 

Static  compression;  823  K 

Samant  et  al.  (1998) 

4.3  GPa 

Full  basal 

Numerical;  single  crystal 

Zhang  et  al.  (2005) 

7.5  GPa 

Full  prism 

Numerical;  single  crystal 

Zhang  et  al.  (2005) 

0.014  GPa 

Partial  basal 

Theory;  Tc=WSFlb 

Maeda  et  al.  (1988) 
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Table  6.  Dislocation  properties  (unless  noted  otherwise:  basal  dislocations,  -room  temperature),  (continued) 


Parameter 

Value 

Defect  type 

Remarks 

Reference 

Tp 

1.24  GPa 

Full  basal  60° 

Theory;  Peierls-Nabarro  stress 

Nabarro  (1947) 

0.060  GPa 

Partial  30° 

tp  =  2  K  exp 

'  -2k  Kd  ^ 

;  — ~1 

1  A 

0.014  GPa 

Partial  90° 

l  &  ) 

0.70  GPa 

Full  prism  edge 

m 

3.0±0.7 

Partial  basal 

Rate  sensitivity  of  yield  stress 

Samant  et  al.  (1998) 

Q 

2.5  eV 

Partial  basal 

Activation  energy;  null  bias 

Samant  et  al.  (1998) 

0.27  eV 

Partial  basal 

Activation  energy;  y=0.5A/cnr 

Galeckas  et  al.  (2002) 

Table  7.  Point  defect  properties. 


Parameter 

Value 

Type 

Remarks 

Reference 

Ef 

3.4  eV 

C 

Vacancy  formation  energy 

Bernstein  et  al. 

(2005) 

5.3  eV 

Si 

sn 

-0.0066  nnT 

C 

Relaxation  volume  (infinite  medium) 

Eshelby  (1954) 

-0.0082  nm3 

Si 

8E2  =  -[3Q^EFl{2iu)]in 
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10° 
10-1  - 
10-2  - 
10-3  - 
10-4h 
j?  10-5h 

o 

0  10-6 1- 


I 

< 

0 

CD 


O 

> 


10-7 

10-8b 

10-9 

10- 


Full  basal  edge  B'=4 
Full  basal  screw  //  =  3.4 
Full  basal  edge  B'  =2.5 
Prism  loop  edge  B'  =4 
Prism  loop  edge  S' =3.4 
Partial  basal  edge  6'  =2.5 
Partial  basal  screw  //'  =  3.4 
Partial  basal  edge  6' =4 


[Legend  and  curves  in  matching 
sequence  of  descending  value 

of  volume  change] 


109  101 


1012  1013  101 


101 


1016  1017  101 


Dislocation  density  [m  ] 


Vacancy  concentration  [ppm] 


Figure  1.  Predicted  volume  changes  from  dislocations  and  vacancies  in  a-SiC. 
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4.  Conclusions 


A  non-linear  continuum  theory  has  been  developed  for  modeling  the  coupled  electromechanical 
behavior  of  anisotropic  elastic-plastic  ceramic  single  crystals  with  defects.  Governing  equations 
and  thermodynamic  definitions  of  physical  properties  have  been  derived  and  concisely  presented. 
Comprehensive  tables  of  compiled  (from  values  in  the  literature)  and  computed  physical 
properties  of  6//-SiC  have  been  presented  to  serve  as  reference  material  to  support  future  meso- 
scale  modeling  efforts  of  ceramic  material  systems. 
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